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ABSTFACT 

Aircraft and Landaat data were uaed with 
coon>uter-alded techniques to del.neate soils 
patterns of a field of 40 ha In a transition zone 
between soils developed under deciduous forest 
and those developed under prairie vegetation. 

Two coeputer-aldad classification techniques , 
aupervlsed and nonsupervlsed, were employed In 
classifying soils of the study area. The means 
and covarlanca matrix statlstlca were obtained 
for every cluster or soil class through the sta- 
tistics algorithm. Each cluster of aircraft and 
Landsat data was Identified and assigned to a 
specific aoll type by correlating the cluster 
soil patterns with a standard aolla map of the 
test site which was prepared as a part of the 
ground observation task. A sampling grid plan 
was usad tt select a training .#et for a super- 
vised classification of the aircraft MSS data. 

The spectral soil patterns revealed In the clas- 
sifications from aircraft and satellite MSS data 
resembled the general patterns of the soils of 
the conventionally prepared soil map. The spa- 
tial resolution of the aircraft scanner was ade- 
quate to recognize each soil type boundary In 
the test site. However, the limited spatial 
resolution of the satellite scanner made It 
difficult to delineate those soil features with 
widths less r'.uin the spatlel resolution of the 
scanner. On the contrary those soil patterns 
which were broad enough to exceed the spatial 
resolution of the Landsat scanner were delineated 
very well. 


of 1:1,000,000, but areas such as floodplains 
ifhlch were too small wejre mapped using 1:250,000 
enlargetoent prints.' 

The general objective of thla Investigation 
was to evaluate and compare the use of computer- 
implemented analysis of multlspect ,1 data from 
aircraft and Landsat scanners to delineate soils 
patterns of one test area In Tippecanoe County 
In Indiana. 


II. STUDY AREA 


A test area of 40 ha was selected In Tippe- 
canoe County, Indiana, In a transition zone be- 
tween soils developed under deciduous hardwocd 
forests and those developed under prairie grasses. 
The soils are within the region of the Alflsols 
but Include some wet Molllsols. The soils In the 
southern half were developed In glacial till with 
less than 40 cm of silt at the surface; whereas 
the soils of the northern half were developed In 
deeper silts. Ttie topography Is level to sloping. 
The following soils are Included In the test arei : 


Reesvllle silt loam 
Cellna silt loam 
Crosby silt loam 
Brookston slit loam 
Brookston silty clay loam 
Ragsdale silty clay loam 
Toronto silt loam 


Aerie Ochraqualf 
Aqulc Hapludalf 
Aerie Ochraqualf 
Typlc Arglaquoll 
Typlc Arglaquoll 
Typlc Arglaquoll 
Udolllc Ochraqualf 


III. PROCEDURES 


I. INTRODUCTION 

Previous experience In reswee multlapectral 
sensing soli studies Irdlcates that the delinea- 
tion of soil boundaries could have limited 
application. Stoner and Horvath demonstrated 
how cultural practices such as plowing and disc- 
ing may affect the multlapectral response of 
surface soils.'' Rrlstof and Zachary also showed 
some limitation In a field being mapped by multl- 
spectral pattern recognition techniques.' Westln 
and Frazee delineated most of the soil associa- 
tion boundaries very well on Landaat Imagery 
using color composite transparencies at the scale 


Multlapectral aircraft data were collected 
on May 6, 1970 by an airborne scanning spec- 
trometer mounted In the University of Michigan 
aircraft at an altitude of 915 m (spatial reso- 
lution 43m' or 0.0043 LANDSAT-2 data 
(spatial resolution of 4500 m' or 0.45 ha) were 
obtained on April 6, 1975 at an altitude of 
915 km. 

A standard soils map was prepared as a part 
of the ground observation task. A sampling grid 
plan was used to select a training set for * 
supervised classification of the aircraft M. 
data. Ten wavelength bands were used In the 
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coaputcr an«lytli of aircraft data. Thaso war# 
0.40-0.44, 0.46-0.48, 0.50-0.52, 0.52-0.55, 0.55- 
0.58, 0.58-0.62, 0.62-0.66, 0.66-0.72, 0.72-0.80, 
and 0.80-1.00 ■Icroaatara. Four wavalangth banda 
wara usad In tha analyala of Landaat-2 data. 

Thaaa wara 0.50-0.60, 0.60-0.70, 0.70-0.80 and 
0.80-1.10 mlcroMtara. 

Two Mthoda of coaputar-alda'I analyala tach- 
nlquaa wara uaad, l.e., auparvlaad and nonaupar- 
vlaad. Tha auparvlaad waa anployad for tha air- 
craft data only. Tha rafaranca aamplaa wara 
aalactad on tha baala of a convantlonal aoll 
aurvay map (Flgura 1 and 2). 



I 1 Raqidaw «icl IBB Toranlo til 

t~ ~1 Brookflon MCI ■■ Crotby Ml 
BroofctKn Ml OB Calmo Ml 
I i Raatavilla Ml 

Flgura 1. Soil Survay Map 
of Teat Site. 

The nonaupervlaed technique waa uaed In both 
aircraft and Landaat-2 data analyala. Tha entire 
teat area of 40 ha from which acanner data were 
collected by aircraft waa aubjected to nonauper- 
vlaed cluaterlng procedurea to obtain fourteen 
apectral or cluater claaaea ualng ten wavelength 
banda. A double number of cluater claaaea waa 
requeated compared to the aeven aoll typea 
occurring In the atandard aoll map to avoid later 
probable Incorrect claaalf Icatlon by the claaalfy- 
polnta algorithm! Since the aame teat area on the 
Landaar-2 data la repreaented by only 88 data 
polnta, a much larger area for cluaterlng waa 
uaed (100 llnea by 100 colunna) . 

To make correlation of remotely senaed data 
with reference eaaler, the Landaat-2 data were 
geometrically corrected before they were uaed.' 

The Landaat-2 data were grouped into 17 cluatera 
ualng four wavelength banda, in an attempt to 


repreaent every type of ground feature. In analy 
ala of both data aourcea (aircraft and aatalllte) 
every aecond data point from every acan line waa 
grouped Into cluatera data having almllar 
apectral charactarlatlca. T.ie Mana and covar- 
iance matrix atatlatlca were obtained for each 
cluater claaa. 

To enhance obaervatlon and to dlacrlmlnate 
the different aoll boundarlaa more eaally, the 
atatlatlcal data were compreaaed Into a ahortar 
format: magnitude of relative reflected energy 

and V/TR ratio (tha aum of relative reflnctad 
energy In the vlalble portion of the apectnai 
divided by the aum of relative reflected energy 
In the reflective IR). 



Figure 2. Computer Soli Claaalf Icatlun 
Map of Teat Site. Legend: (M) Ragsdale 

silty clay loam, (F) Brookston silty clay 
loam, (0) Toronto silt loam, (C) Celliui silt 
loam, (-) Reesvllle silt loum, (/) Crosby 
s(l: loam, (I) Brookston silty loam, (A) ve- 
getation. 

Corresponding statistics In the fotm of 
magnitude and ratio were assigned to each 
cluster class. Baaed on these statistics the 
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Figure 4. Level 1 computer mep of test site 
from slrcraft ratio data. t.egend: 

- “ high ratio; I - medium rntlo; 

F “ low ratio. 

vlded Into high, medluoi, low , and lowest 
groups (Figures 5 and 6) . Using LAKD- 
SAT-2 data, soils of the same area were separated 
Into more levels than with aircraft data, fn the 
first phase, 13 of the 17 cluster clar.' <?s w>'re 
Identified as bare soil by analysis of .ev< 1 I 
statistics. The statistics from the r •uter ana- 
lysis were used In LARSYS merge and gl 'tint pro- 
c;shors to produce computer lesult nap with high, 
i.'>’dlum, and low soil spectral response ^ Figures 


RESULTS and DISCUSSION 

The soils were well separated from other 
non-soil classes of the study site from aircraft 
and satellite MSS data. The reflectance patterns 
of soils at various wavelengths are considerably 
different from all other suterlal on the ground. 

Comparing a standard soil survey nap of the 
AOha teat site with a computer-aided supervised 
tap of aircraft MSS data (Flgurr 2), one can see 
that the spectral patterns revealed In the compu- 
ter classification nap resrnblet' the general pat- 
terns of the soils of the convent 1 nally prepared 
soil nap (Figure 1). Some small areas of Brooks- 
ton silty clay loan arc laapped as Ragsdale soils 
and vice versa. Light-colored Reeaevllle soils 
are mapped very well. Cellna and Crosby soils 
have the same drainage characteristics and rlml- 
lar surface color as Reeaevllle soils. Toronto 
and Brookston silt loam mapped by computer-aided 
techniques are In good agreemen. with the stan- 
dard soil survey nap. The spatial resolution of 
the aircraft scanner was adequate to recognize 
each soil mapping unit In the test site. 


A hierarchical approach for soil cLmsK (ca- 
tion was used In both aircraft i.nd LANDSAT dst.i 
analyals. The general separation of s>'>11b in four 
levels Is bssed only on spectral InfcrmAt l'>n , ob- 
serving the Mgnltude and ratio betwe-en c<> 'h of 
the soil cluster classes separately. Thr ' 
observed with aircraft scanner data an- ipci trslly 
divided on Level I Into high, medium, and low res- 
ponse soils (Figures 3, 4). Level II lu lubdl- 


Flgure 3. Level 1 computer map of test site 7 and 8). Fifty of the dst.- points fell Into 


from aircraft magnitude data. groups of soil with low magnitude and high spec- 
Legend: - • high siagnltode; 1 “ tral ratio values. Thirty points were of medium 
medium magnitude: F • low magnitude. magnitude, and only one data point had high mag- 
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hierarchy of solla was astabllshad for this 
Invest Igat Ion (Tables 1,2,3 and 4). Cluster 
class areas were merged Into two levels of air- 
craft and three levels of landsat-2 data. Level 
I la composed of three categories for aircraft 
and Landaat-2 data. Level II contains seven 
catcscrlee for alrcreft and ala categorlea fer 
Landsat-2 data. Level III consists of ten soil 
categories of Landsat-2 data only. 


After cluaters were grouped Into desired 
soil categories, the Level I, IT and III statis- 
tics were used as training statlstlca for Input 
In the supervised classification approach. An 
overlaid Interpretation technique was used to 
coaipare soil categories on photo enlargements 
made from computer classification maps. 


Table 1. Hierarchy Baaed on Magnitude Developed for Soil Spectral Investigation of Aircraft Data. 



Level I 

Level II 

Response 

NS-Clasa 

Symbols 

Coda 



High A 

323.60 

NS-1/14 

♦ 

Ha 


High 

301.12 

278.65 

NS-2/14 

• 

Hb 


279.87 

High B 

256.04 

NS-3/14 

• 

He 



258.77 

261.50 

NS-4/14 

m 

Hd 

y 


Medium A 

218.55 

NS-5/14 

/ 

Ma 

«J 

Medium 

219.32 

220.10 

NS-6/14 

I 

Mb 

9C 

219.32 

Medium B 

198.35 

NS-7/14 

j 

Me 

> 


202.40 

206.46 

NS-8/14 

z 

Hd 

S 


Low A 

191.06 

NS-9/14 

c 

La 



184.35 

17/. 65 

NS-10/14 

0 

Lb 


Low 

Low B 

164.52 

NS-12/14 

A 

Lc 


173.99 

156.03 

147.54 

NS-13/14 

H 

Ld 



Low C 
131.15 

131.15 

NS-14/14 

F 

Le 

Table 2. 

Hierarchy Based on 

Ratios Developed 

1 for Soil 

Spectral Investigations of Aircraft 

Data. 


Level I 

Level II 

Resprnse 

NS-Class 

Symbols 

Code 


High 

1.38 

High A 
1.46 

1.46 

NS-4/14 

+ 

Ha 


1.33 

NS-1/14 

. 

Hb 


High B 

1.38 

NS-5/14 

- 

He 



1.36 

1.37 

NS-7/14 

- 

Hd 



Medium A 

1.30 

NS-2/14 

/ 

Ma 



1.29 

1.29 

NS-3/14 

I 

Mb 


Medium 






«» 

1.27 


1.26 

NS-6/14 

J 

Me 



Medium B 

1.26 

NS-9/14 

Z 

Md 

*T 

e 


1.26 

1.27 

NS-10/14 

C 

Me 



Low A 

1.24 

NS-8/14 

0 

La 



1.24 

1.25 

NS-12/14 

A 

Lb 


Low 

Low B 

1.21 

NS-13/14 

H 

Lc 


1.18 

1.21 







Low C 
1.03 

1.03 

NS-14/14 

r 

Ld 
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Table 3. Hierarchy Baaad on Magnitude Davalopad for Soli Sr>actral Invaatlgatlona of Landaat>2 Data 



Uval t 

Uvel II 

Level III 

IU.eponae 

HS-Claaa Synbola 

Code 



High A 

High A 

153.03 

HS-1/17 

■f 

Ha 


High 

143.96 

143.96 

137.72 

HS-3/17 

e 

He 


132.86 

High B 

High B 

122.92 

NS-5/17 


Hb 



122.92 

122.92 








Medlun C 

133.66 

NS-7/17 

• 

Me 



Hedlua A 

113.46 







113.46 

Medlun D 

U3.13 

hS-8/17 

/ 

Md 

•p 

Madlu* 


113.13 





m 

W 

108.30 


Medlun E 

lPi.02 

NS-9/17 

I 

Me 

m 


Medliw B 

105.02 





• 






1 


104.47 

Medlun F 
103.44 

103.44 

HS-11/17 

J 

Mf 




Low C 

97.11 

NS-14/17 

z 

U 



Low A 

97.11 







96.57 

Uw D 

96.03 

NS-12/17 

c 

U 




96.03 






Low 








88.38 


Low E 

89.80 

NS-15/17 

0 

Le 




86.06 

86.73 

NS-13/17 

A 

Lg 



Low B 


82.14 

NS-16/17 

H 

Lh 



83.80 

Uw F 

71.30 

NS-17/17 

F 

Lf 




71.30 






Table 4. Hierarchy Baaed 

on Ratloa Developed 

for Soil 

Spectral Invaatlgatlona of Undaat-2 Data. 


Uvel 1 

Uvel II 

Level 

Ill Reaponee NS-riaaa 

Synbola 

Code 


* 

High A 

High A 

High A 

1.74* NS-17/17 

• 

Ha 


High 

(1.52) 

1.74 

1./4 

High B 

1.51 NS-14/17 

/ 

Hb 



High B 

High B 

High C 

1.44 NS-14/17 

• 

He 



1.45 

1.45 

High D 

1.40 NS-16/17 

+ 

Hd 




Medlun C 
1.35 

Medlun C 

1.35 NS-7/17 

L 

Me 



Medlun A 

Medlun D 

Medlun D 

1.33 NS-9/17 


Md 

• 

u 


1.33 

1.33 





p 

4i 

« 

Medlun 


Medlun E 

Medlun E 

1.31 NS-5/17 

0 

Me 

M 

(1.30) 


laBl 





1 



Medlun F 

Medlun F 

1.27 NS-3/17 

J 

Mf 

0 


Medlun B 

1.28 







1.26 

Medlun G 
1.24 

Medlun C 

1. 1 NS-13/17 

I 

Mg 




Low C 

Low C 

1.20 NS-1/17 

8 

Lc 



Low A 

1.20 






Low 

(1.16) 

1.19 

Low D 
1.19 

Low D 

1.19 NS-12/17 

A 

Ld 



Low B 

Uw B 

Uw E 

1.15 NS-8/17 

4 

U 



1.13 

1.13 

Uw F 

1.12 NS-11/17 

F 

Lf 
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Figure 5. Level 2 computer map of test eite 

from aircraft magnitude data. Leg- 
end: • ■ high magnitude; - « 

medium magnitude; 0 • low magnitude; 
4 ■ lower magnitude; F “ lowest 
magnitude. 



Figure 6. Level 2 computer map of test site 
from aircraft ratio data. Legend: 

• ” high ratio; I “ medium ratio; 

0 ■ low ratio; 4 • lotrer ratio; F « 
lowest ratio. 


nltude value. Figures 9 and 10 show Level II 
classifications with five soil sub-groups In which 
separation Is based on spectral response In the 
form of magnitude mi ratio between visible and 
reflective IR portions of the spectrum. Again, 
the majority of the data points went Into two 
classes of low, and leas Into medium and high 
reflective soils. A more detailed spectral dif- 
ferentiation of the soils Is obtained In Level III 
where tne soil tast area la broken down Into 
eleven spectral groups (Figures 11 and 12). 

In order to achieve greater spectral con- 
trast, the study area was extended and more clus- 
ter classes were Introduced In the analysis of 
the Landsat data. This procedure contributed very 
little In separating the two low reflective soils, 
namely Ragsdale and Brookston silty clay loaios. 
This may be expected, because the Brookston soil 
series consists of very poorly drained, nearly 


level soils with a very dark gray surface, while 
the Ragsdale soil series consists of deep, dark- 
colored, poorly drained soils with a black silty 
clay loam surface layer. To obtain a better 
separation of the soil series, 21 samples of dark- 
colored soils were evaluated with two data points 
per sample and 16 samples of medium and light- 
colored were evaluated with two data points each. 
The samples were sorted Into an array from lowest 
to highest spectral response levels. The samples 
were grouped Into eight spectral classes based on 
magnitude of reflectance and reflectance ratio. 
These classes were used as reference classes In 
machine-aided classification. The automcted LAND- 
SAT classification map (Figure 13) was compared with 
an aircraft classification map. This comparison 
revealed that the large and homogeneous areas of 
soils could be delineated from LANDSAT-2 data. 
Small mapping areas are merged together in 
larger LANDSAT classification areas, or they 
are added to areas with similar spectral proper- 
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Flcuie /. Level 1 coaputer mep of test site 
froB Lendeet magnitude data. Leg- 
end; M - farmatead; - • high mag- 
nitude; 0 ■ medium magnitude; F • 
lov magnitude. 


tlea. The apatial resolution of the aatelllta 
scanner system la such that la la not adequate 
for delineation ot soil mapping units with an 
extension of only a few hectares. 
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000 — 

□ -00 

□ “0 F 

00 — F 

0 — 000 — 

Q-030 

- 0 - 000 " 

-00000-F 


Figure 8. Level 1 computer map of test site 
from Landaat ratio data. Legend; 
M ■ farmstead; - - high ratio; 0 ' 
medium ratio; F - low ratio. 
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Flgura 9. Level 2 conputer lup of teat alte 
fron Landaat nagnltuda data. Lag- 
end: M ■ farnatead; / • high mag- 

nitude; I “ medium auignltude; 0 ■ 
lov magnitude; A “ lower magnitude; 
F “ lowaat Mgnltude. 
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Figure 10. Laval 2 computer map of teat alte 
from Landaat ratio data. Legend: 

N • farmstead; - « high ratio; / ■ 
medium ratio; 0 “ low ratio; 4 ■ 
lower ratio; F ■ lowest ratio. 


1977 A^achne Processing o/t RemoteiK Sensed Dcrta Symposium 


I 



M/AOO=OZ 
= 10001*0 
lOOZOIOC 
OAAI=IAH 
I AIIOAFA 
»00»0HAJ 
000 »» 00 J 
IOAI-=OA 
=011 1000 
0*01 =IOA 
01=1 1 lOZ 


Plgur* 11 • L«v«l 3 conputer map of teat alte 
from Landaat magnitude deta. Leg- 
end: N = farmateud; - * highest 

magnitude; / ” high A magnitude; 

• ■ high B magnitude; I ■ medium A 
magnitude; J “ medium B magnitude; 

0 - low A magnitude; Z = low B mag- 
nitude; A “ low C magnitude; C “ 
low D magiiltude; H ■ low R magni- 
tude; P “ low F magnitude. 
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Figure 12. Level 3 costputer map of test site 
from Landsat ratio data. Legend: 

H - farmstead; - - highest A ratio; 
+ - highest B ratio; • ■ high A 
ratio; / ■ high B ratio; L = medium 
A ratio; Z ■ medium B ratio; 0 ■ 
medium C ratio; I ■ medium E ratio; 
A • low B ratio; 4 • low C ratio; 

F ■ low D ratio. 
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Figure 13. Co«i>utar up of test alt* fro« 

LANDSAT data related to atandard 
aoll up. Legend! ■•/,!• 
Reeaeville allt loan; 8>Croaby 
allt loan; O.F^Ragadale allty 
clay loam; A-Brookaton allty 
clay loaai; Z>Brookaton allt 
loan. 


ferance between two nolle depended on the depth of 
allty aurfaca utarlal. 

Autyala of LANDSAT data with co«puter-alded 
technlquaa ahowed that It waa not ponalble to find 
apectrally ho«c -/*r»''ua aoll featurea of the aavan 
aoll nerlea on the AOha teat alte on the digital 
dlaplay or on a plcture*prlnt up. On the other 
hand, cluatarlng technlquaa could be uaed on an 
extended teat area to group apectrally alallar 
data polnta Into cluater clannea. Cluatar claaa 
ntatlatlca In the form of ugnltude and ratio 
nerve aa a baala for grouping. The level claaaea 
arc then related to the aoll patternn. In ac«e 
caaen the LANDSAT HSS data ware not adequate for 
raaolvlng aoll featurea with wldtha lean than that 
of the acanner aynr..i'a npatlal renolutlon (ap- 
proxlutely 70a). Thona aoll patterna which were 
broad enough to exceed the apatlal renolutlon of 
the LANDSAT acanner were delineated very well by 
apectral analynla. 

Typically, the total field of view Increaaea 
aa the altitude of the data collection ayaten In- 
creaaea. However, luge renolutlon decreaaea aa 
altitude Increaaea, ao there la lean detail avail- 
able fro* high altltudea. Advantagea and dlaad- 
vantagea of both high and low reaolutlon acanner 
ayateu muat be taken Into account If coaiputk.T- 
aldad analyala technlquaa are to be uaed aa a 
baala for aoll curvey. 
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CONCLUSIONS 

Computer-aided analyala technlquaa uaed with 
aircraft MSS data aho%red that the apatlal reao- 
lutlon was sufficient to recognize each noil 
mapping unit of the rest site. Seme difficulties 
occurred where different soil nerlea were Intri- 
cately mixed and this mixture showed aa a sepa- 
rate spectral mapping unit, or where the dlf- 
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